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Abstract: We consider interface fluctuations on a two-dimensional layered lattice 
where the couplings follow a hierarchical sequence. This problem is equivalent to the 
diffusion process of a quantum particle in the presence of a one-dimensional hierarchical 
potential. According to a modifled Harris criterion this type of perturbation is relevant 
and one expects anomalous fluctuating behavior. By transfer-matrix techniques and by 
an exact renormalization group transformation we have obtained analytical results for the 
interface fluctuation exponents, which are discontinuous at the homogeneous lattice limit. 
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I. Introduction 



Recently there is a growing interest in such natural and artificial systems which are or- 
ganized in a hierarchical way. Examples can be found in economical organizations [1] and 
stock-market exchanges [2], in geological processes before major earthquakes [3], in stud- 
ies of relaxation phenomena of proteins [4], spin glasses [5] and computer architectures [6]. 
Theoretically much effort have been devoted to understand the linear dynamics (i.e. the 
diffusion process) in a system with hierarchically organized energy barriers. According 
to numerical[7,8] and exact[9,10] results the diffusion in such systems can be anomalous 
(which is often called as " ultradiffusion" [7] ) , furthermore in several models there is a dy- 
namical phase transition [8] separating regions with normal and anomalous diffusion. For 
a comprehensive review on the subject see Ref[ll]. 

Another subject of theoretical interest is the properties of (static) phase transitions 
on hierarchical lattices. For these and other non-periodic (quasi-periodic or more gen- 
erally aperiodic) systems a relevance-irrelevance criterion has recently been proposed[12] 
on the analogy of the Harris criterion[13] for random magnets. The cross-over exponent 
corresponding to a non-periodic perturbation is given by: 

^ = l + uD{n-l) (1) 

in terms of the v correlation length exponent of the unperturbed system and the wandering 
exponent of the sequence 0[14]. Here D denotes the number of coordinates on which the 
couplings depend, c.f. for a layered system D = 1. The perturbation is then expected to be 
relevant (irrelevant) if $ > ($ < 0), which was indeed found in a series of exact studies 
on two-dimensional layered Ising models[15,16]. For marginal sequences, where $ = 0, 
continuously varying critical exponents and anisotropic scaling behavior was observed[17]. 

As far as the critical behavior on hierarchical lattices is concerned mainly the two- 
dimensional layered Ising model with a one-dimensional Huberman-Kerszberg (HK) se- 
quence[7] and the corresponding Ising quantum chain were studied. In numerical[18] and 
exact [19,20] calculations non- universal critical behavior were found in accordance with the 
vanishing cross-over exponent in eq(l), which follows from the fact that the fluctuation 
exponent of the HK sequence is f2 = 0[20]. 

In this paper we consider the interface fluctuation problem on a layered lattice, where 
the couplings between the layers follow the HK hierarchical sequence. As far as interface 
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wandering on non-periodic lattices are concerned the work by Henley and Lipowsky[21] has 
to be mentioned, who considered the interface roughening in two-dimensional quasicrystals. 
On a layered lattice with Fibonacci-type quasi-periodicity non-universal interface fluctu- 
ations were observed, with a continuously varying interface wandering exponent. This 
behavior is again in accord with the relevance-irrelevance criterion, since with f2 = — 1 
and V = vji_ = 1/2 the cross-over exponent in eq(l) is $ = 0. In our problem, on the HK 
lattice = 0, thus $ = 1/2 > and the perturbation is relevant. Therefore one expects 
anomalous interface fluctuations on this lattice. 

The structure of the paper is the following. We deflne the model in Sec. II. The 
results of the transfer matrix calculations and that of an exact renormalization group 
(RG) transformation are presented in Sections III. and IV., respectively. The results are 
discussed in the flnal section. 

II. Formalism 

We consider a diagonally layered ferromagnetic spin model (c.f. the Ising model) on 
the square lattice with hierarchically organized interactions. The couplings in the /i-th 
diagonal — Jh/ksT are selected from a set {kq, ki, K2, ■■■) and Kn — tikq, such that 

Kh = Kn , /i = 2^(2m+l). (2) 

This type of structure of the couplings (Fig.l), which shows the typical features of ultra- 
metric topology [5] was introduced by Huberman and Kerszberg[7] following the work in 
Ref[l]. 

The boundary spins on the (1,1) surfaces are flxed in different orientations (Fig.l) and 
we are interested in the fluctuations of the interface separating the (+) and (-) regions. 
The interface is considered as a continuous structurless string and complicated interface 
configurations, such as overhangs and bubbles are omitted. It is generally accepted that to 
study interfacial fluctuations it is sufficient to keep only Solid-on-Solid (SOS) type interface 
configurations. In this so called SOS model the interface is geometrically represented by a 
directed walk or polymer [22]. 

In the SOS model the interface is characterised by its h{x) height at site x and the 
interfacial energy is specified by the Hamiltonian: 

-H/kBT^Y.'^Ky,^-) ' (3) 

X 
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where surface effects are omitted. The thermodynamic properties of the interface are 
conveniently studied in the transfer matrix formahsm[23,24]. For our model the transfer 
matrix in the x-direction, parallel with the boundaries, is given by: 

Th,i = 5h,i-ie-''''^+5h,i+ie-^''' . (4) 

Here according to eq(2) the matrix-elements are from a set (eo, ei, €2, ...) and the ratio 
of successive terms is constant: en+i/cn = -R < 1. For the homogeneous system = 1, 
whereas for hierarchical lattices R measures the strength of inhomogeneity. The interface is 
not likely to visit sites with a matrix-element e„, n ^ 1, since the corresponding probability 
is weighted by a factor of i?". 

The interfacial free energy cr and the longitudinal correlation length ^y, which is 
measured parallel with the boundaries, are given in terms of the leading and the next-to 
leading eigenvalues of the transfer matrix Aq and Ai as: 

cr=-logAo (5) 

and 

e|7^ = log(Ao/Ai) . (6) 
The fluctuations of the interface grow on a power law scale: 

< [/i(0) - /i(a;)]' >~ , (7) 

where w is the wandering or fluctuation exponent, which is = 1/2 for homogeneous 
two-dimensional systems [22]. 

Another quantity of interest is the probability Pq{x) that the interface after a;-steps 
has the same position, i.e. h{Qi) = h{x). For a walk or diffusion problem, where x plays the 
role of the time, Pq{x) is the autocorrelation function, which has the asymptotic behavior 
Pq{x) ~ x~'^ . For homogeneous two-dimensional lattices 7 = 1/2 and generally w = 7(8]. 
It could be shown by slightly modifying the derivation in Ref. [8] that the autocorrelation 
function averaged over the starting positions of the interface can be expressed through the 
spectrum of the transfer matrix as: 
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where g{X) = i/L'^-6{X — Aj) is the density of states and L denotes the width of the 
system in the h direction, thus it is the dimension of the transfer matrix. 

The eigenvalues of the transfer matrix are dense at the top of the spectrum and one 
can develop a scaling theory in terms of these critical eigenvalues. We consider a critical 
level Aj of a system with a finite width L and denote by AAj = Aq — Aj its difference from 
the top of the spectrum. Changing lengths by a factor of b=2, i.e. with L' = L/2 the z-th 
eigenvalue will be A^, and the difference (AA^)' will scale with a factor of b^^, thus 

{AXiy = 2y^AXi , (9) 

where yx is the gap exponent. We stress that the statement in eq(9), that all critical 
levels scale with the same factor is a scaling hypothesis, which will be verified by actual 
calculations in the following sections. 

Using eq(9) the transformation law for the density of states is given by 

^(AA) = 2y^-'g'[{AXy] (10) 

which is compatible with a power law dependence of the density of states at the top of the 
spectrum: 

^(AA) ~ (AA)^/^^-^ . (11) 

Now putting this expression into eq(8) and evaluating the autocorrelation function one 
gets 7 = 1/yA. 

From the scaling behavior of the spectrum in eq(9) one obtains for the finite size 
corrections to the largest eigenvalues: 

Ao - A,(L) ~ L-y^ , (12) 

thus from eqs(6) and (12) the longitudinal correlation length is ~ L^^. In a finite system 
the correlation length perpendicular to the (1,1) surface is limited by the width of the strip 
L, therefore the interface wandering exponent in eq(7), which can be alternatively 
defined as ~ is given by: 

w = l/yx . (13) 

Thus indeed = 7, as expected from scaling considerations. 

In the following we calculate the interface fluctuations on the HK lattice by two 
methods. First we study numerically the spectrum of the transfer matrix, verify the 
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validity of the scaling hypothesis and determine the interfacial tension and the wandering 
exponent. Then we apply an exact renormalization group transformation and calculate 
analytical expressions for the critical exponents. 

III. Numerical Study of the Transfer Matrix 

The transfer matrix of the interface problem in eq(4) is tridiagonal and could be diag- 
onalised by powerful methods [25]. In the specific problem, however, due to the hierarchical 
structure of the transfer matrix one can implement a very fast algorithm to calculate the 
roots of the corresponding determinant. 

We consider a finite system of size L = 2' and express the corresponding determinant 
D(2') by two subdeterminants of sizes 2^"^ and 2^~^ — 1, respectively, in the form: 

D{2^) = D(2'-i)D(2'-i) - D(2'-i - l)D(2'-i - l)eti • (14a) 

The symmetric determinant D{2^ — 2) of size 2' — 2, which is obtained from Z)(2') by 
leaving out the first and last rows and columns, can be similarly expressed as: 

D{2^ -2) = D{2^-^ -1)D{2^-^ -1) - D{2^-^ -2)D{2^-^ -2)ef_:^ . (146) 

Finally: 

D{2^ -1) ^ D{2^-^)D{2^-^ -1) - D{2^-^ -1)D{2^-^ -2)e^_-^ . (14c) 

These relations supplemented with D{1) = -A, D{2) = - eg and D{2) = A^ - 
define a fast procedure to calculate the value of the determinant for very large sizes. For 
example we could treat with this method slightly perturbed systems with i? pa 1 up to 
sizes L = 2^0 - 2^°. 

The largest eigenvalues calculated by this method all have the same type of finite 
size dependence, thus the scaling hypothesis in Sec. II is indeed satisfied. The leading 
eigenvalues calculated on the largest finite lattices are accurate at least up to 10-12 digits. 
The gap exponents, describing the finite size dependence of Xi{L) in eq(12), however could 
be obtained from the raw data with a comperatively smaller accuracy, up to 5 digits. In 
this case to increase accuracy we used sequence extrapolation methods, such as the van 
den Broeck and Schwartz and the Bulirsch and Stoer methods[26]. 
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The leading eigenvalue of the transfer matrix, which is connected to the interfacial 
tension in eq(5) and the extrapolated values of the the interface wandering exponent are 
listed on Table 1. One can see that both the leading eigenvalue and the wandering exponent 
are monotonically decreasing as R goes from one to zero. In the limit — > the interfacial 
tension in eq(5) together with the wandering exponent go to zero, which is due to the 
fact that the system tends to be separated into disconnected parts. More interesting is 
the behavior of the wandering exponent around the homogeneous lattice point. As the 
value of R is lowered below one the wandering exponent jumps by a finite amount of 
Aw = 0.0432799 from w = 1/2. In the renormalization group language such type of 
behavior corresponds to a relevant perturbation, which brings the system into another 
stable fixed point. In the next section we shall explicitly construct the RG transformation 
and determine exactly the wandering exponent. 

IV. Renormalization Group Calculation 

We are going to study the scaling behavior of the largest eigenvalues of the transfer 
matrix in eq(4), which satisfy the second order difference equation: 



where in the thermodynamic limit the boundary terms are omitted. The structure of the 
couplings which are connected to in eq(4) are shown on Fig.l. To construct an exact 

recursion we decimate out those sites, which are connected to a ni coupling or equivalently 
to an ei matrix-element (denoted by crosses on Fig.l). We note, that the same type of 
decimation was used by Maritan and Stella in their study of the diffusion problem on the 
HK lattice[10]. One can see that after a decimation step the (eo,ei,eo) triplet will play 
the role of the renormalized eg, whereas the other couplings will renormalize as = e^+i, 
keeping the value of R and together the structure of the transfer matrix unchanged. 

Performing the RG transformation first we denote the two neighbouring sites to be 
decimated out by i and i + 1 and express and V'i+i ^is: 



(15) 



(16) 
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where A = eoA/(A^ — ef) and B = eoei/(A^ — ef). Then the difference equations in terms 
of the remaining, non-decimated spins have the same form as in eq(15), provided the 
eigenvalue and the couphngs transforms as: 

w _ A - Aeo 

< = ^ , n=l,2,... (17) 

and Cq = cq. Thus the ratio of the sequence remains invariant R' = R, as expected. As a 
consequence in the RG transformation besides A only one coupling, say ei is enough to be 
considered and the RG transformation can be written as a two-parameter recursion: 

X2 -.2 ^2 

eoei 

e[ = R^^ , (18) 
Co 

where cq is the input value of the largest matrix-element. 

The physically relevant fixed point of the transformation with A > is given by: 



R f XV Vl-R + R'^ 



6q J 1 — R \,6q J 1 — R 



(19) 



which is stable for < < 1. The eigenvalues of the linearised fixed point transformation 
are roots of a quadratic equation and they are given as: 



Ai2 = — + -R+-± 
^'^ R 2 



1 r. 

— + R+- 

R 2 



1/2 

(20) 



The leading eigenvalue Ai > 1 determines the scaling behavior of the spectrum of the 
transfer matrix and the y\ scaling dimension is given by: 

The second eigenvalue of the RG transformation is A2 < 1 and the corresponding scaling 
field is irrelevant, thus the fixed point in eq(21) is attractive and governs the critical 
properties of the physical model with ei = Reo. It is seen from eq(19) that the fixed point 
of anomalous interface fluctuations does not exist at the homogeneous point i? = 1, where 
the fluctuations are characterised by the normal wandering exponent w = 1/2. Comparing 
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the analytical results for = 1 /y\ with those obtained by finite size calculations in Table 
1, we can say that the numerical results are indeed very accurate they correspond to that 
in eq(21) at least up to six digits. 

V. Discussion 



In this paper we studied interface fluctuations on a layered hierarchical lattice. The 
perturbation caused by inhomogeneous couplings is relevant according to a linear stability 
analysis and the observed interface fluctuations are indeed anomalous. The wandering 
exponent is a monotonically decreasing function of R and discontinuous at i? = 1. The 
fact that w{R) < 1/2 can be understood, since the interface preferentially stays on kq 
lines and the probability to visit a Kn line is rapidly decreasing with n. Consequently the 
interface fluctuations are damped by the inhomogeneously distributed couplings. 

One can estimate w{B) in the limit i? — > 0, when the probability of a large interface 
fluctuation of height = 2" is primary given by Pn ~ e^, i.e. by the probability to have 
one step on the Kn line. For such a fluctuation the interface approximately takes x ~ 

~ steps, thus the wandering exponent in leading order is: w{R) = — log 2/ log i?, 
which corresponds to the aymptotic behavior of the analytical result in eq(22). We note 
that in the ^ limit the interface fluctuations can be described by a Markovian process 
and then our problem is equivalent to the diffusion of a particle in a hierarchical lattice, 
as studied in Ref[7-ll]. 

The HK-sequence used in this paper can be generalized, by having a general v-djcy 
character [27], instead of z/ = 2 used in eq(2). Then one has in eq(2) h = R'^{um + /j,), with 
fj, = 1,2, ...,1/ — 1. According to our numerical and analytical investigations for v = 3 and 
u = A the main characteristics of interface fluctuations remain the same as for = 2: the 
wandering exponent has a jump at i? = 1 and varies with R. For v = 3 we obtained the 
analytical result: 

log 3 



^u=z = 2 ( 4 + i? + 1 1 + 



log Ai.=3 

^ X 2 1 1/2 



^ ■ (22) 

As mentioned before the problem studied in this paper is related to the diffusion pro- 
cess on hierarchical lattices [11]. Our problem, however, can be formulated as the quantum- 
mechanical diffusion process of a particle which is represented by a wave packet and placed 
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on a one-dimensional HK-potential. Then x and h{x) correspond to the t time and the po- 
sition of the particle at the given time step, respectively, while the transfer matrix describes 
time evaluation. According to our results in a one-dimensional hierarchical potential the 
width of the wave packet will grow in time anomalously as t'^^^\ 

Our final remark concerns some similarities of our results to that of interface fluctua- 
tions in a repulsive, inhomogeneous surface potential, decaying as ~ where I measures 
the distance from the surface[28]. In two-dimensions for uj < 2 the perturbation is relevant 
and the interface wandering exponent takes the anomalous value: w = 1/u > 1/2 [29]. 
In this problem, however, the perturbation is confined to the surface, furthermore the 
wandering exponent is continuous at a; = 2. 
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R Ao/eo w = 1/yx 

1. 2. 0.5 

0.999 1.99800894 0.4567199 

0.9 1.82853274 0.4551092 

0.75 1.62218648 0.4451438 

0.5 1.35286081 0.4004540 

0.25 1.14948652 0.3110577 

0.1 1.05381456 0.2272971 

0.001 1.00050038 0.0911867 



Table 1 



The leading eigenvalue and the corresponding interface fluctuation exponent from 
numerical diagonalization of the transfer matrix for different values of the hierarchical 
parameter. 
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Figure captions 



Fig.l: Structurless interface on a diagonally layered square lattice. The values of the cou- 
plings, which follow the hierarchical HK sequence in eq(2) are indicated below. Sites 
to be decimated out in the RG transformation are marked by X. 
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